EPJ manuscript No. 

(will be inserted by the editor) 



O 



o Fluctuation spectra of free and supported membrane pairs 



(N 
O 



OO 



o 

CO 

-4-* 



o 
o 



> 
m 
in 

O 
+-> 



i 

C 

o 



x 

S3 



Rolf-Jiirgen Merath 1 ' 2 and Udo Seifert 3 

Max-Planck-Institut fur Metallforschung, Heisenbergstrafie 3, 70569 Stuttgart, Germany 

2 Institut fur Theoretische und Angewandte Physik, Universitat Stuttgart, Pfaffenwaldring 57, 70569 Stuttgart, Germany 

3 II. Institut fur Theoretische Physik, Universitat Stuttgart, Pfaffenwaldring 57, 70550 Stuttgart, Germany 

February 5, 2008 

Abstract. Fluctuation spectra of fluid compound membrane systems are calculated. The systems addressed 
contain two (or more) almost parallel membranes that are connected by harmonic tethers or by a continu- 
ous, harmonic confining potential. Additionally, such a compound system can be attached to a supporting 
substrate. We compare quasi-analytical results for tethers with analytical results for corresponding contin- 
uous models and investigate under what circumstances the discrete nature of the tethers actually influences 
the fluctuations. A tethered, supported membrane pair with similar bending rigidities and stiff tethers can 
possess a nonmonotonic fluctuation spectrum with a maximum. A nonmonotonic spectrum with a maxi- 
mum and a minimum can occur for an either free or supported membrane pair of rather different bending 
rigidities and for stiff tethers. Typical membrane displacements are calculated for supported membrane 
pairs with discrete or continuous interacting potentials. Thereby an estimate of how close the constituent 
two membranes and the substrate typically approach each other is given. For a supported membrane pair 
with discrete or continuous interactions, the typical displacements of each membrane are altered with re- 
spect to a single supported membrane, where those of the membrane near the substrate are diminished 
and those of the membrane further away are enhanced. 

PACS. 87.16. Dg Membranes, bilayers, and vesicles - 82.70.-y Disperse systems, complex fluids 



1 Introduction 

Research on biological membranes has benefitted consider- 
ably from studying well-defined model systems pQ. Mem- 
branes attached to a substrate can be investigated with 
surface-sensitive techniques (e.g., reflection interference 
contrast microscopy, fluorescence microscopy, atomic force 
microscopy, surface plasmon resonance, and quartz crystal 
microbalance) . However, for solid-supported membranes, 
the substrate tends to demobilize and denaturalize mem- 
brane proteins. Hence, systems with a biocompatible sur- 
face are intended to mimick the biological membrane in 
its natural environment [2]. A promising and meanwhile 
widely used approach is to tether the membrane to the 
substrate, ensuring a larger distance between the mem- 
brane and the solid while keeping the membrane confined. 
Recent experimental achievements and developments per- 
mit (i) to micro- and nano-structure substrates in order to 
configure the pattern of attachment sites [3j[4j[5] , (ii) to tai- 
lor the length and stiffness of the tether molecules (or solid 
pillars) 6,7,905], and (iii) to adapt the linker chemistry 
to different substrates [10l[TT][T2] . Different structures of 
composite membranes are generally accompanied by dif- 
ferent thermal fluctuation spectra, which have previously 
been measured for red blood cells [13] and for a bilayer 
membrane interacting with a laterally homogeneous sub- 
strate [Til. 



Theoretically thermal shape undulations of mem- 
branes have been addressed for various model systems. 
Fluctuation spectra of supported membranes have been 
calculated for a single membrane connected to a fixed sup- 
port. This connection was established either by a conti- 
nuous harmonic confinement [15] or by tethers at discrete 
attachment sites [16lfT7lfT8lH9] . For the latter case of local- 
ized interactions, nonmonotonic fluctuation spectra have 
been observed. The substrate structure has been consid- 
ered in the context of adhering supported membranes |20j . 
The attachment of a membrane to an elastic meshwork 
has also been studied [2T]. Nethertheless, so far, fluctua- 
tions of an underlying support have not yet been included 
quantitatively into the analysis of membrane fluctuations. 

In a next step towards more biocompatible [22] soft 
interfaces and to more realistic biomembrane models one 
could attach a membrane not directly to a substrate but 
rather to a second, supported membrane. Since presently 
several experimental efforts are going into this direction 
[23ll24ll25] . we expect supported membrane systems with 
more than one fluctuating membrane to be constructed in 
the near future (e.g., by utilizing vesicles tethered to a 
supported membrane [26] and subsequent vesicle fusion). 
These systems could be used (i) to study diffusion within 
a membrane in a native-like environment, (ii) to tune, to a 
certain degree, the fluctuation properties of the constitu- 
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Fig. 1. Side-view sketches of the model systems addressed in this article. These compound membrane systems are built up by 
almost parallel fluid membranes in an aqueous solution, potentially situated on a solid support (substrate). Temperature-driven 
shape fluctuations are influenced by the interactions between the membranes and/or the substrate. The models A-C comprise a 
continuous harmonic confining potential between the constituent membranes and — if present — the substrate, whereas discrete 
tethers constitute the connections in the models D-F. In the models A and D a single membrane is attached at a substrate, in 
B and E two membranes are linked together but are mobile in the vertical direction as a whole, and in C and F a membrane 
pair is additionally attached to a substrate. The single membrane models A and D have already been discussed in our previous 
work [18] and provide a basis for comparison of the different models in limiting cases. 



ent membranes in order to investigate fluctuation-depen- 
dent diffusion of membrane proteins [171127] or to suppress 
membrane fluctuations in order to avoid undesired contact 
[2"B] , and (iii) to prolongate the lifetime of biosensors [55] . 
Contemporary passivations of the substrate usually last 
only a few days and the lower membrane could then act as 
an alternative passivation that prevents the upper mem- 
brane to contact the substrate. 

The purpose of this article is twofold: (i) The fluc- 
tuation behavior (in Fourier space and in real space) of a 
tethered two-membrane system on a substrate is analyzed. 
The discreteness of the tether sites demands numerical 
calculations even though the result can be expressed in a 
quasi-analytical way. (ii) Smearing out the discreteness of 
the tethers results in a continuous confining potential. For 
this system the fluctuation spectrum for arbitrary ratios 
of membrane bending rigidities as well as the standard 
deviation in real-space for equal bending rigidities are dis- 
played analytically. Typical membrane displacements as a 
function of the ratio of bending rigidities are shown. Con- 
tinuously interacting membrane systems are interesting on 
its own, nethertheless we present it here as a limiting case 
to which the results of the discrete model are compared. 
The systems considered here are supposed to be far from 
unbinding transitions 30,31], thus steric repulsion can be 
subsumed in an effective harmonic binding potential 32, 

Esg. 

This article is organized as follows: In section [2] we 
introduce the models addressed in our work. In section 
[3] the models containing continuous harmonic interac- 
tions are considered, and fluctuation spectra are derived 
analytically. In section 2] we approach systems compris- 
ing membranes connected by discrete tethers and derive 
quasi-analytic solutions for the corresponding fluctuation 
spectrum and root-mean-square displacement (RMSD). 
The results from the continuous and discrete model sys- 
tems are compared. A summarizing perspective is given 



in section \5\ Finally, the case of three membranes with 
continuous interactions is considered in the Appendix. 



2 Models 

Several model systems of fluid compound membranes are 
analyzed in this article. In Fig.[TJ an overview of the differ- 
ent cases is given. The diversity of possible configurations 
is classified with respect to three structural features: We 
distinguish between (i) discrete tethers or continuous con- 
finement, (ii) one or two (or even more) membranes and 
(iii) "free" (i. e., mobile as a whole) or supported systems. 

The behavior of a particular model can be understood 
qualitatively by considering limiting cases that correspond 
to other models from the above scheme. An interesting 
question is, under what circumstances the discrete nature 
of the tethers determines the fluctuation behavior signif- 
icantly and makes numerical effort necessary. Thus, we 
compare the class of model systems with discrete teth- 
ers to corresponding systems with a continuous harmonic 
potential acting on elongation differences of the involved 
membranes. 

Notations are explained exemplarily for model F: The 
basal plane is given by a rectangle with side lengths L x 
and L yi and the bending rigidities of the two membranes 
are denoted K\ and K2- The two equally large membrane 
surfaces are parametrized by h (1) (r) = h {1) (x.y) (upper 
membrane) and h (2) (r) = h (2) (x,y) (lower membrane), 
where h {1) = and h (2) = are the rest positions of the 
upper and lower membrane, respectively. In the rest posi- 
tion the two membranes are flat, parallel and separated by 
the rest length l of the N equally long, vertical tethers. 
Each tether between the two membranes is described by 
its lateral site r a and its spring constant K a . Each site 
of the substrate-membrane tethers is given by 

(s) 

and its spring constant by K a ■ In contrast to the dis- 
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Crete tethers, continuous harmonic interactions involve an 
elastic parameter 7 that can be related to a comparable 
spring constant (cf. equation ([50)) ). 

In our theoretical model sytems the membranes can, in 
principle, penetrate each other and the substrate. Hence 
the results are only applicable to those systems for which 
such fluctuations are irrelevant. This is the case if the rest 
length lo of the tethers is considerably larger than typical 
membrane elongations which can be estimated a posteriori 
by the root-mean-square displacement. 



3 Membranes with continuous interacting 
potential 

We first study systems with interactions between a mem- 
brane and a nearly parallel object like another membrane, 
a substrate or an underlying cytoskeleton. Here, these 
interactions will be approximated by a continuous, har- 
monic, confining potential. Such models with continuous 
interactions also serve as comparison with the discrete 
tether-linked models addressed in the next section. 



3.1 Model B: Membrane pair with continuous 
interaction 

As the simplest two-membrane model system the case of 
two membranes connected by a continuous binding po- 
tential (model B, cf. Fig. Q} is studied. The Hamiltonian 
of the system consists of a Helfrich bending energy term 
for each membrane and the potential energy due to the 
interaction, 

H m = [ L dx / i dy{^i[V 2 fc (1) (r)] 2 + — [V 2 /i (2) (r)] 2 
Jo Jo 2 2 

+ l[h^(r)-h^(r)} 2 } . (1) 

The index 1 or 2 refers to the upper or lower membrane, 
respectively. The parameter 7 determines the strength of 
the harmonic potential. A Fourier ansatz for both heigth 
profiles, 

h»\r) = ^ k ^e ik ' r , (2) 

with i G {1, 2}, k = (k x ,k y ) and k x . y = (0, ±1, ±2, . . .) • 
(2ir/L Xy y), leads to 

H^ = ±L x L y J2 k [(«i|k| 4 + 7 )l^ 1) | 2 

+ ( K2 |k| 4 + 7 )|/ l < 2, | 2 -2 7 Re(/ 1 < 1) C*)], (3) 

where Re denotes the real part. Separating the Fourier 
coefficients into real and imaginary parts, 



results in 



2 L * L vI2 k 



[( Kl |k| 4 + 7 )( a L 1,2 + &: 
+ Mk| 4 + 7 )(aif + 



(1)2 
k 



-2 7 KX>+6<X 



(5) 



Since the height profiles are real, the Fourier coefficients 
fulfill the conditions a ( \ — and 6^' k = — £>£\ so that 
we chose the coefficients belonging to the wave vectors {q} 
and as independent variables. In order to avoid divergent 
translational modes, we exclude the translational mode 
dp of the lower membrane. Or, in other words, the verti- 
cal position a { Q ' serves as frame of reference. This choice 
is motivated by related model systems containing a sub- 
strate: Since a substrate on the one hand calls for acting 
as frame of reference and on the other hand corresponds to 
an infinitely rigid membrane, we do not turn to the center 
of mass as frame of reference but to the spatial average of 
the lower membrane height profile. Thus the Hamiltonian 
as a functional of the independent variables is 



H<- 



x± ^y 



[( Kl |q| 4 + 7 )(< )2 + 6< 1)2 ) 
( K2 |q| 4 + 7 )(< )2 + &< 2 > 2 ) 

2 7 (<<+«) 



\ L x L yl a o )2 - (6) 



Integrating out the degrees of freedom of the lower 
membrane, {a^'} and {&q'}, leads to an effective Hamil- 
tonian for the upper membrane, 



(7) 



0H ™ (8) 



Applying the auxiliary formula 

00 „ 

dz (2) e -/3L x L B ( K2 |q| 4 + 7 ) z ( 2 )-/3L,L H (-2 7 ) z ( 1 > 2 ( 2 > 

■00 

P U|q| 4 + 7 7 ) 
for each integral in equation (JU) gives 



(9) 



(10) 



After a suitable shift of the energy zero-level, the effective 
Hamiltonian reads 



-^cff = L x L y ^Ki|q| 4 + 7 



7 



"k 



for i g {1,2} 



(4) 



q 



K 2 |q| 4 +7 

T T (!) 



(11) 
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The fluctuation spectrum of the upper membrane can 
now be read off. For non- vanishing wave vectors it is given 

by 



k^Ol 



*k / t- \ « k / - 2 ( a 
ksT / L x L y 



) 



Kl |k| 4 + 7 - 
kBT/L x Ly 



T 



K 2 |k| 4 + 7 



K X k 4 +7eff(A:) 
with an effective potential parameter 



7off(fc) 



1 

7 



h: 2 



k 4 



(12) 
(13) 

(14) 



and the abbreviation k = |k|. This formula does not apply 
for the translational mode of the upper membrane, i.e., 
for the spectral origin k = 0, since Oq 2) was set to zero 
manually. Instead, at the origin the spectrum attains the 
value 



<id 2 > = <ar> - 



kBT/L x L v 

7 



(15) 



If the translational mode clq 



of the lower membrane 
was permitted (i. e., if it was regarded as a degree of free- 
dom), the spectrum would be described correctly for all 
wave vectors by 



K\ 2 : 



kBT/L x Ly 
ki k 4 + 7cff(fc) 



(16) 



featuring a divergence at the origin. 

Examining "f e g from equation (|14|) a striking analogy 
is emerging. j c g depends on 7 and ^k 4 like the effective 
spring constant of two springs in series depends on the con- 
stituent spring constants. Thus, the lower membrane af- 
fects the upper membrane like a second continuous spring, 
arranged in series to the connecting continuous springs. In 
contrast to the "real" continuous spring parameter 7, the 
strength K^k 4 of the "virtual" continuous springs depends 
on the wavelength. 

The limiting cases for 7 e ff are 



7eff 



7 for 7 <C K2 k 4 
K2 k 4 for 7 > K2 k 4 , 



(17) 



and thus for the spectrum 



ksT/LxLy 



K\ 2 \ 



for k > ( 7 /k 2 ) 1/4 



K\ k 4 + 7 

^ k « (7/.2) 1 / 4 



(18) 



Consequently, for large wave vectors, the case of a single 
membrane linked to a substrate by continuous springs is 
recovered, while for small wave vectors, the upper mem- 
brane behaves as if two free membranes were glued to- 
gether exhibiting an effective bending rigidity K\ + «2. 



x 1 



(1) 2 membranes, K1/K2 — 1 

(2) limiting case: 1/[(1+ 

(3) limiting case: + 1) 

(4) free membrane, 1/zf 

(5) 2 membranes, — 10 

(6) 2 membranes, /C1/K2 — 10 -1 

(7) 2 membranes, /C1/K2 — 10 -2 

(8) ^ limiting case: 1/[(1+ lfl 2 ): 

(9) 2 membranes, K1//C2 — 10 - " 1 



Fig. 2. Fluctuation spectra for model B. A membrane is con- 
tinuously, harmonically coupled to an almost parallel second 
(lower) membrane. For several values of the ratio ki / K2 of the 
bending rigidities, the scaled spectrum is plotted versus the re- 
duced wave vector z\ containing the bending rigidity k\ of the 
upper membrane and the stiffness parameter 7 of the harmonic 
potential. At the origin the scaled spectrum always diverges, 
if the translational mode of the lower membrane is permitted, 
whereas it equals 1, if this mode is artificially excluded or the 
lower membrane is regarded as the center of the frame of ref- 
erence. The limiting curves for large or small z\, respectively, 
intersect at z\ = (jci/k^) 1 / 4 , indicating the position of the 
cross-over region. Evidently, the spectrum becomes similar to 
the spectrum of a single membrane attached to a substrate for 
K2 — > 00, and the spectrum of a free membrane is recovered 
for K2 — » 0. 



Introducing the reduced wave vector 
zx = k£i , 
with the persistence length 

6 = (kV 7 ) 1/4 , 

the spectrum can be expressed as 



K\ 2 : 



k B T 



(19) 



(20) 



(21) 



Its limiting cases are 



k B T 
7 L^Ly 



Z? + l 



fo"i»(|) I/4 
fo-i«G|) 1/4 



(22) 



The spectrum of the upper membrane and its limiting 
cases are displayed in Fig. [2] for varying values of ki/k2- 



3.2 Model C: Membrane pair on a substrate, with 
continuous interactions 

If the membrane pair from the previous subsection (model 
B) is attached to a substrate by a continuous harmonic 
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interaction, the Hamiltonian reads 





H^= / dx / dy{^[V 2 /i (1 '(r)] 2 + ^[V 2 /i (2 '(r)] 2 
^ (2) (r)] 2 + |[^(r)] 2 } . (23) 



This situation is equivalent to a system of three mem- 
branes with continuous interacting potentials (cf. the App- 
endix), in which the bending rigidity of the lowest mem- 
brane is set to infinity, K3 — > 00, and the translation mode 
Oq is excluded. The fluctuation spectra of the upper mem- 
brane. 



(l^ 1> | 2 > = 

Ki k 4 + 71 

and of the lower membrane, 



k B T/L x L y 



7i 



(24) 



«2 k 4 + 71 + 72 



(K } \ 2 ) 



kBT/L x Ly 



K.-2 



k 4 



(25) 



72 



7i 



follow from equations (HHJ) and (|.14p . respectively, and also 
hold at the origin k = 0. 

The spectra above can be written as 



I ft' 1 " 2 



n 



k B T 



1 



and 



id 2 : 



L xLy 71 zf + l- 



k B T 

L X Ly 72 



(26) 



7i 



4 + 1 + 



(27) 



"KT"~4 

K 2 2 



ITT 

72 



(28) 



(29) 

is stretched to- 



with the reduced wave vectors (i £ {1, 2}) 

Zi = k d 
and the persistence lengths 

6 = ( Kl / 7i ) 1/4 • 

On increasing ji, the spectrum ( \h { k \* 
wards higher wave vectors and its amplitude is reduced. 
In contrast to the tether models discussed below, the con- 
tinuous harmonic potential does not involve any physical 
length, and the length scale of the system is merely given 
by the persistence lengths £j. For the simplified case of 
equally strong interacting potentials, i.e., 7 = 71 = 72, 
these two types of spectra are displayed for various values 
of ki/k 2 in Fig. [3] 

The real-space fluctuations of this system are consider- 
ed now. In general, the mean-square displacement (MSD) 
of the membrane i results from the correlations of the 
Fourier coefficients, 



[^(r)] 2 )=E(^^)e i(k+k,) 



(30) 



(a) 2 



(1) Kl/«2 00 («2 - 

(2) K!/k 2 = 2 

(3) Kilm = 1 

(4) Kl /K2 = 0.5 

(5) Ki/m = ltT 1 

(6) ki/k 2 = 10- 2 

(7) Ki/m = 10- 1 

(8) Ki/m = 10- 5 

(9) Ki/k 2 -> (k 2 — 



— r~ 
0.5 



(b) 1 




(1) Kl/te^O 

(2) Ki/K2 = 10- 

(3) Ki/m = 0.5 

(4) Kxlm = 1 

(5) Ki/rc2 = 2 

(6) Ki/m = 10 

(7) K1/K2 = 10 2 

(8) ki/k 2 = 10 3 

(9) Kl /iC2 = 10 5 

(10) ki/k 2 



Fig. 3. Fluctuation spectra of the upper (a) and lower (b) 
membrane in model C, i. e., in a system consisting of two mem- 
branes coupled among each other and to a substrate via a con- 
tinuous interacting potential. The scaled spectra (cf. equations 
(|26|) and (J2TJ) ) are plotted versus the reduced wave vector z\ 
and Z2 (cf. equation (|28[) h respectively. Both elastic constants 
are set equal, 71 = 72 • 

(a) The spectrum of the upper membrane recovers the spec- 
trum of a single membrane attached to a substrate with 7 = 
I/271 for ft 2 — > and with 7 = 71 for ft 2 — > 00 (except at the 
origin, since the translational mode is not affected by the 
bending rigidity ft2 and obeys ( \h^ 2 ) = 2 x ksT / L x L y -y\). — 

(b) The lower membrane complies with a single membrane at a 
substrate with 7 = 71 for fti — > and with 7 = 2 71 for fti — > 00 
(except at the origin, since fti does not hinder the translational 
mode from maintaining ( |fto 2) | 2 ) — 1 x ksT / LxLy^)- 



Since the present system with continuous interacting po- 
tential is translationally invariant, it is 



(W) = (l^l 2 >^- k 



(31) 



Consequently, the MSD of the upper (i = 1) and lower 
(i = 2) membrane can be expressed as 



ia£T: 



(32) 



k,k' 



We consider the special case of equally strong interact- 
ing potentials, i.e., 7 = 71 = 72. After simple manipula- 
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3/%/5 



1 

2/^5 



1.414 
1.342 



1 

0.894 



-- 0.707 



10" 



10" 



io- 



1 



10- 



10 4 



10 6 



Fig. 4. Scaling functions b m (n\j ' K2) and 6 <2 '(«;i/k2) of the up- 
per or lower membrane, respectively, in model C. They describe 
the amplification of the MSD in a system of two membranes 
with a continuous binding potential among each other and to- 
wards a substrate with respect to a single membrane attached 
to a substrate. Here, the two involved continuous potentials 
are assumed to be equally strong. 



tions, the MSD of the upper membrane reads 



([/^(r)] 2 ) = J^= x - / ,1,- 



2 £L 

K2 



VT^T tt Jo z 4 + (1 + 2^).* 2 + t 
k B T 



(33) 



where the first factor corresponds to the MSD of one single 
membrane linked to a substrate by an interacting potential 
with parameter 7 and the scaling function describes 
the influence of the lower on the upper membrane. One 
can read off the two limiting cases, 



& (1 Vi/^ 2 ) 



1 for k 2 -> 00 
V2 for k 2 -► . {M) 



Analogously, the MSD of the lower membrane reads 




x b (2 \ Kl /K 2 ) 



(2 + ff)- 2 + f 
(35) 



The two limiting cases of the scaling function 6 <2) are 





1 



for Ki 
for K\ — > 00 . 



(36) 



For the case of equally rigid membranes, k\ = k 2 , the 
scaling functions take the values = 3/v5 ~ 1.342 

and 6 <2) (1) = 2/\/5 ~ 0.894. The scaling functions and 
are shown in Fig.[4j Compared to the single-membrane 
case, the fluctuations of the upper membrane are increased 
and those of the lower membrane are decreased. This re- 
sult has an impact on related experimental systems, since 
the risk of undesired adsorption of the membrane at the 
substrate is reduced by introducing another membrane on 
top. 



4 Two membranes linked by (discrete) tethers 

After having analyzed the models in which the membranes 
are subject to a continuous harmonic confinement, we 
address the models containing localized connections via 
tethers. First in subsection I4.il the basic formulae are 
derived for model E (cf. Fig. [I}, where a coupled mem- 
brane pair can move freely in the vertical direction. Sub- 
sequently a solid support, composed of tethers as well, 
is added (model F) and incorporated into the calculation 
scheme fsubsection 14. 2j) . Then the results are discussed 
with respect to fluctuation spectra (subsection I4.3[) and 
typical membrane elongations (subsection [ 



4.1 Model E: Basic calculation 

For the system of two membranes linked by tethers (model 
E, cf. Fig.[T|) the Hamiltonian consists of a bending energy 
term for each membrane and the potential energy of the 
N tethers due to membrane elongations at the tether sites 

ff« ee J L dxJ% [V 2 ^'(r)] 2 + ^[V 2 ^(r)] 2 } 
+ E A \ ^ (I) (r a )-^ (2) (r Q )f ■ (37) 

^ ^ OL— 1 Z 

The involved quantities are explained in section^ In or- 
der to avoid a divergence of the translational modes, we 
preliminarily exclude the translational mode of the lower 
membrane or, in other words, switch over to the frame of 
reverence given by the spatial mean height of the lower 
membrane (cf. subsection 13. 1|L 

Inserting Fourier expansions for both membrane pro- 
files, 

fc«(r) ee £ft£V k -' (38) 

k 

with i S {1, 2}, k= {k x , k y ) and k x>y = (0, ±1, ±2, . . .) • 
(2ir/L x , y ), and translating to independent variables [17] . 
the Hamiltonian can be expressed in the following form 

m, 

H = d T Md ee £d t M t>t ,d t , . (39) 
t,t> 

The vector d of independent system variables is built up 
by the vectors c (i) of independent variables of both mem- 
branes, 

'c (1) \ 

jm) ■ (4°) 

Their components c ( r ' ] (with r = 0, 1, 2, . . .) can be grouped 
into three sectors, 



c<" T ee (§fcg>, {Reft™}, {Im/i«} 



(41) 



where q runs through all independent, non- vanishing wave 
vectors and Re and Im denote the real and imaginary 
part, respectively. Note that, if the translation mode of a 
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membrane is manually excluded (to avoid its divergence), 
the corresponding vector and matrix elements have to be 
discarded. 

The symmetric, positive definite matrix M consists of 
four quadrants and can be written as 



M = 



/M (1 > H 
\ H M (; 



(42) 



if the set of possible wave vectors are the same for both 
membranes. The submatrices M (1) and M (2) represent a 
single membrane tethered to a substrate [18] , whereas the 
matrix H contains the interaction of Fourier modes of both 
membranes. These matrices can be simplified to 



and 



with 



and 



N 



r< • m a r , (43) 



N 



to" ■ m" 



(44) 



Q = l 



E { *> = 





&i L x L y |q(r)| 4 for r > 



for r = , 

4 _n ( 4 5) 



with n™ defined analogously to m" (cf. equation (|4"6")l ) 

(s) 

employing K a instead of K a , i. e., n" 



r 'Ka—tK 



(b). 



If the solid support is built up by a continuous interact- 
ing potential instead (cf. the last term in equation (j2"3"l) ~). 
the associated parameter 7 must be added to E {2) , i.e., 
Ef> -» £< 2 > + 7 . 



4.3 Fluctuation spectra 

Quasi-analytical results for two membranes (discretely) 
tethered between each other and possibly also to a sub- 
strate are presented and compared to the correspond- 
ing models with continuous interactions. First fluctuation 
spectra for the different cases are analyzed, then real-space 
elongations are examined. 

In order to make the continuous harmonic potential 
comparable to an array of equally strong discrete tethers, 
we replace the stiffness parameter 7 by the parameters of 
the discrete model, that lead to the same effective spring 
constant for planar membrane elongations, 



1 T ^ NK K 

L x L y ^-'0=1 a L x L y A 2 



(50) 



Then the reduced wave vector z, used in the continuous 
cases, is related to kA via 



^/2~K~ a for r = (1 st sector) 

///, = <; yJ2K a cos[q(r) • r Q ] for the 2 nd sector (46) 
y/2K a sin[q(r) • r a ] for the 3 rd sector . 

The inverse matrix M _1 provides the correlations of 
all Fourier coefficients, 



kA 



(d t d t 



(47) 



thus in particular the fluctuation spectrum | 2 ) of 

each membrane can be extracted. 



4.2 Model F: Additional solid support 

The extension to the membrane system additionally teth- 
ered to a solid support (model F, cf. Fig. [T]) is straight- 
forward. Let ./V"( s ) vertical, harmonic tethers of strength 

Ka at lateral sites r„ attach the lower membrane to the 
substrate. The corresponding Hamiltonian reads 



^— 'a=l 



(s) 



(48) 



The additional term finds its corresponcance in the cal- 
culation scheme by a modification of the submatrix M <2) 
(cf. equation 



N 



M; 2 ;, = E ( r 2> ■ 5 rr , + > m° ■ m a r , + > < • r# (49) 



a=l 



a=l 



(KA 2 /n) 1 / i 



(51) 



Except for the free membrane case (and pathological 
cases in the tether models), the fluctuation spectrum ex- 
hibits a relative extremum at the spectral origin. We do 
not count this trivial extremum and, in the following, con- 
sider extrema only for non- vanishing wave vectors. 



4.3.1 Fluctuation spectra of model E (without substrate) 

Fluctuation spectra of a tethered membrane pair (model 
E, cf. Fig. [T]) are considered. This compound membrane is 
mobile as a whole since there is no contact to a substrate. 
As a choice for the tethers connecting the two almost par- 
allel membrane sheets, we consider a quadratic array of 
equally strong tethers, K a = K. The translational mode 
of the lower membrane is excluded (cf. section l4~Tj) , and we 
focus on the upper membrane and its spectrum ( l^k'l 2 )■ 
Several spectra belonging to model E are shown and com- 
pared to limiting cases in Fig. [5] 

Tethered membrane spectra are anisotropic for large 
enough spring constants (cf. Fig. 2 in Ref. [TB]), while 
spectra associated to continuous interactions are isotropic. 
Since the spectrum for a membrane system mobile as a 
whole covers several orders of magnitude in the interest- 
ing region, azimuthal averages are vulnerable to numerical 
artefacts. Thus the spectra in Fig. [5] are evaluated along 
a main axis in Fourier space (e.g., k = (fc^O)), which is 
related to a main axis (r = (x,0)) of the quadratic tether 
meshwork in real-space. 
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We distinguish between four different, partially over- 
lapping regions in the parameter space constructed by the 
two parameters K1//C2 and KA/k%: (i) a soft lower mem- 
brane, (ii) a rigid lower membrane, (iii) weak tethers, and 
(iv) relatively stiff tethers and a stiff lower membrane. 

The relationship between the different membrane mod- 
els (cf. Fig. Q} reveals the following two trivial limit- 
ing cases. The case (i) of a soft lower membrane (i.e., 
K1//C2 - > °o) corresponds to a free membrane (except at 
the origin, due to choice of the frame of reference). Con- 
versely, for (ii) a rigid lower membrane (i.e., K1/K2 — » 0) 
the behavior of a single membrane attached to a substrate 
is recovered. 

(iii) For weak tethers (i.e., for KA 2 /ki < f0 2 ) the 
spectrum (l^'l 2 ) is well approximated by the formula 
for two membranes jointed by a continuous potential (cf. 
section l3Tj) . It is reasonable that the spectrum of a mem- 
brane tethered to a second membrane lies above the cor- 
responding spectrum for attachment to a substrate (that 
complies with K2 — > 00). Thus it is plausible that the 
continuous description cannot be valid for wave vectors 
for which the "continuous spectrum" lies below the sin- 
gle membrane "pinning" spectrum. The intersection of 
the continuous and the pinning spectrum marks a cross- 
over region between a small wave vector region, where the 
continuous formula is valid, and a large wave vector re- 
gion, where the spectrum is well described by the pinning 
spectrum. Thus for long-wave undulations the discrete- 
ness of the tether meshwork becomes less important, and 
for short-wave undulations the lower membrane acts like a 
substrate. Note that for K — > 00 the continuous spectrum 
is equivalent for all wave vectors (cf. equation (|18p) to the 
behavior of a single membrane with an effective bending 
rigidity K ett = «i + n 2 . 

The case (iv) of relatively stiff tethers and a stiff lower 
membrane comprises a nonmonotonic (two-dimensional) 
spectrum ( \ h^ | 2 ) with a relative minimum and a relative 
maximum. The monotonic and the nonmonotonic regime 
are separated by a line in (ki/«2, K A 2 / Ki)-sp&ce connect- 
ing two points, firstly {k\/k2 = c\ ~ ICP 2 , KA 2 /ki — > 00) 
and secondly {ki/ks — > 00, KA 2 / k\ — c-i — I0 2 ). The 
nonmonotonic region is located on the side of this "coex- 
istence line" which possesses smaller K1/K2 or/and larger 
KA 2 /k\ values. Hence, in other words, nonmonotonicity 
in model E occurs only for relatively stiff tethers and two 
membranes possessing very different bending rigidities. 

All spectra for model E share the qualitative beha- 
vior at the spectral origin (k = 0) where the transla- 
tional mode of the lower membrane has been excluded to 
avoid a divergence at the origin. Within numerical accu- 
racy we can state that ( |/io'| 2 ) for discrete tethers coin- 
cides with the spectrum of a single membrane tethered 
to a substrate. Thus, all discrete spectra in parts (a) and 
(b), respectively, of Fig. [5] take the same value at the ori- 
gin and then jump to a high value (which is related to 
the smallest wave vector larger than zero and therefore to 
the extensions of the basal plane). This can be understood 
by a simple consideration: For an equidistant meshwork of 
equally stiff tethers let any Fourier mode (except the a pri- 
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Fig. 5. Fluctuation spectra for model E, i.e., of a membrane 
attached to a second (lower) quasi-parallel membrane by a 
quadratic meshwork of equally stiff tethers. The frame of refer- 
ence is the spatial average of the lower membrane heigth profile, 
thus all tethered membrane spectra in (a) and (b), respectively, 
coincide at the origin. Cuts through the spectra along a main 
axis are displayed. 

(a) Rigid tethers (K — > 00). The spectrum (|/ik'| 2 ) exhibits 
different regimes for varying values of Ki/«2. A single mem- 
brane "pinning" spectrum is found for k\/k2 — > 0. For increas- 
ing ratio the relative minimum moves away from the 
origin and is less pronounced. The spectrum remains nonmono- 
tonic up to K1/K2 — 1CF 2 (curve 2), where a transition between 
nonmonotonicity and a monotonicity occurs. A cross-over re- 
gion between the (long wavelength) behavior as if the discrete 
tethers were a continuous potential and the (short wavelength) 
behavior as if the lower membrane was completely stiff (like a 
substrate) is situated around the intersection of both limit- 
ing curves. The expression for the continuous interaction spec- 
trum takes a simple form in the present case, namely that of 
a free membrane with a modified bending rigidity. Finally, for 
K1/K2 — > 00 the upper membrane fluctuates like a free mem- 
brane. — (b) Relatively strong tethers. In the displayed case of 
KA 2 /ki — 10 2 the (two-dimensional) single membrane spec- 
trum is at the transition between monotonic and nonmonotonic 
behavior. In comparison to the situation in (a) , here the conti- 
nuous interaction spectrum represents a better approximation 
for wave vectors close to the cross-over region. 



ori forbidden translational mode) of the lower membrane 
push the upper membrane away at some lateral position 
(x, y) = {x,y). Then there is always a symmetric posi- 
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tion (x, y) = [L x — x, L y — y), where a tether pulls in the 
other direction. Thus the net effect on the lateral spatial 
average (that equals Zip ) of the upper membrane height 
profile vanishes, i.e., undulations of the lower membrane 
do not influence ( l^o'l 2 )■ 



4.3.2 Fluctuation spectra of model F (with substrate) 

After the freely moving two-membrane system, we now in- 
vestigate the model with additional solid support tethered 
to the lower membrane (model F, cf. Fig.[T]). For simplic- 
ity we assume the additional substrate-membrane tether 
meshwork to be identical to the quadratic membrane- 
membrane tether meshwork, i.e., = N, — K, 



and ri s) = r a . 

In the tethered membrane model the tether distance 
A imposes a natural length scale. In contrast to the con- 
tinuous interaction model (cf. equation (|19|) ). the spectra 
cannot be expressed in terms of a reduced wave vector 
that incorporates the spring constant. The spectrum of 
the upper or lower membrane can be expressed in terms 
of a scaling function with i £ {1,2}, 



(K\ 2 ) 



kBT xfi) (n 1 

KiN \ K 2 



KA 2 



,k x A,k y A) . (52) 



We proceed as follows. For both the upper and the 
lower membrane we distinguish three cases depending on 
the parameters ki/k 2 and KA 2 / 'm: (a) rigid tethers (i. e., 
K — > oo), (b) a vanishing bending rigidity of the other 
(i.e., not focussed) membrane, and (c) similar bending 
rigidities. Case (c) itself contains three regimes, (i) stiff 
tethers, (ii) intermediate spring constants, and (iii) weak 
tethers. Only the upper membrane exhibits an additional 
regime (d) for stiff tethers and a stiff lower membrane, 
where a rich nonmonotonic spectrum is found. 

(a) For rigid tethers (K — > oo) the spectrum of a 
single membrane pinned to a substrate is recovered inde- 
pendently from the two bending rigidities due to the two 
tether meshworks situated directly on top of each other 
(in "phase"). 

(b) For a small bending rigidity n 2 of the lower mem- 
brane (i.e., large values of n\/n 2 ) bending modes of the 
lower membrane are energetically favored. Therefore the 
upper membrane behaves like a single membrane attached 
at a substrate by discrete tethers with an effective spring 
constant K eft = K/2, since from the point of view of the 
upper membrane the two tether meshworks are aligned in 
series. Analogously for small k.\/k 2 the lower membrane 
behaves like a single membrane attached to a substrate by 
discrete tethers with spring constant K, since the lower 
membrane feels the two tether meshworks like arranged 
in parallel. 

(c) For intermediate values of ki/k 2 (i.e., ki/k 2 ~ 
0(1) ) the spectra of both membranes can be grouped into 
three regimes according to the tether stiffness, (i) For large 
spring constants (KA 2 /Ki > 10 3 ) each spectrum coincides 
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Fig. 6. Fluctuation spectra of the upper (a) and lower (b) 
membrane in model F, i. e., for two identical membranes (/ti = 
K2) tethered among each other and to a substrate. Both tether 
meshworks form quadratic arrays of equally stiff tethers situ- 
ated directly on top of each other. The two-dimensional spec- 
tra are averaged with respect to the azimuthal angle in the 
(k x , fc H )-plane and plotted versus the radial part k. 
(a) For stiff tethers the spectrum of the upper membrane be- 
haves like a single membrane spectrum with an effective spring 
constant K e!{ = K/2 (curves 1 and 2). It converges to the spec- 
trum of a two membrane model with continuous interaction for 
small spring constants (curves 7 and 8). In between, i. e., in the 
regime of intermediate values of the spring constant (curve 4) , 
slight deviations from those limiting cases (curves 5 and 6) 
appear for undulations modes with wavelengths around a few 
tether distances. — (b) The lower membrane features the same 
spring constant regimes as the upper membrane, but the limit- 
ing case of a single membrane (i. e., for stiff tethers) possesses 
the same spring constant (K ctt = K) as the original tethers, 
since for the lower membrane the two tether meshworks are 
arranged in parallel and not in series. 



with the discrete single membrane spectrum with an ef- 
fective spring constant K c{{ . Focussing on the upper mem- 
brane necessitates K cH = K/2, while for the lower mem- 
brane K en = K applies. Consequently, a nonmonotonic 
spectrum of the upper or lower membrane, respectively, is 
observed for K ott A 2 / Ki > 10 2 [TB]- (ii) For an intermedi- 
ate tether stiffness small, suppressive deviations from the 
corresponding single membrane spectrum appear around 
kA ~ 7T. Finally, for (iii) weak tethers (KA 2 /^ < 10) 
the formulae for a supported membrane pair with conti- 
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kA 

Fig. 7. Nonmonotonic fluctuation spectra for unlike bending 
rigidities in model F. Two relative extrema (apart from the ori- 
gin) appear in the spectrum of the upper membrane in model 
F (cf. Fig. [1} , if the lower membrane and the tethers are rel- 
atively stiff. The new relative minimum for small wave vec- 
tors has evolved from the relative minimum at the origin of 
a nonmonotonic single membrane spectrum, since long wave- 
length undulations of the lower membrane enlarge the upper 
membrane's spectrum around the origin. Azimuthal averages 
of the spectrum of the upper membrane are shown for rela- 
tively stiff tethers (KA 2 /ki — 10 3 ) and various bending rigid- 
ity ratios Ki/k2- For comparison, also related spectra of sin- 
gle tethered membranes (curves 1,2, and 8) are displayed. For 
K1/K2 — > the spectrum resembles that of a single tethered 
membrane with the same spring constant (except at the ori- 
gin). For K1/H2 — » oo the spectrum of a single tethered mem- 
brane with half the spring constant is recovered. For any value 
of Ki/«2 the spectrum is situated between these two limiting 
cases, but always maintains the same value at the origin. The 
structure at kA ~ n is an artefact from the azimuthal av- 
eraging and does not appear in the original, two-dimensional 
spectra. 



nuous interactions (cf. subsection l3.2p hold. Spectra of the 
upper and lower membrane for K1/K2 = 1 for the above 
mentionned three spring constant regimes are displayed 
and compared to the limiting cases in Fig. [6] 

(d) A qualitatively new effect arises for stiff tethers 
(at the order of the critical spring constant KA 2 /n\ cAt ~ 
100, at which a single tethered membrane possesses a 
nonmonotonic spectrum) and for a relatively stiff lower 
membrane: Then the upper membrane behaves like a sin- 
gle tethered membrane (with spring constant K and not 
K/2) only for intermediate and large wave vectors. More- 
over, for small wave vectors the — despite of the large 
K2 — active translational mode of the lower membrane 
affects the upper membrane and thus increases the spec- 
trum around the origin. For a sufficiently stiff lower mem- 
brane (ki/k 2 < 10~ 3 ) the spectrum of the upper mem- 
brane exhibits a relative minimum at a small wave vector 
in addition to the familiar relative maximum. This new 
type of nonmonotonicity is shown in Fig. [7l In the case 
of a single tethered membrane with stiff tethers a relative 
minimum is localized at the origin, but it moves to a non- 
vanishing wave vector for two membranes and suitable 



parameter values. In the limiting case of K1/K2 — > oo the 
spectrum of the upper membrane looks like a single teth- 
ered membrane spectrum except at the origin where the 
spectrum jumps to a higher value. Of course, very strong 
tethers hide this effect and yield the well-known pinning 
spectrum. 

After the detailed discussion of the various regimes we 
qualitatively illustrate how the two membranes influence 
each other and lead to deviations from the behavior of 
a single tethered membrane. The influence of the lower 
on the upper membrane can be understood as follows: 
Undulations of the lower membrane displace the attach- 
ment sites of the upper tethers. With increasing bend- 
ing elasticity of the lower membrane the probability of 
these attachment site fluctuations decreases. The concept 
of an effective spring constant only works well for negligi- 
ble suppression of such displacements, so that two tethers 
on top of each other can be regarded as two springs in se- 
ries. Since the effective (in general wave- vector-dependent) 
spring constant is smaller than K (and equal to K/2 in the 
limiting case of k% — > 0) , the spectrum of the upper mem- 
brane always exceeds the corresponding single membrane 
spectrum. 

On the other hand, the influence of the upper on the 
lower membrane is slightly different: Here the limiting case 
for vanishing bending rigidity of the other membrane cor- 
responds to a single tethered membrane spectrum with 
an effective spring constant equal to K. From the point of 
view of the lower membrane the two tether meshworks are 
arranged in parallel. A non- vanishing bending rigidity n\ 
leads to an effective (in general wave-vector-dependent) 
spring constant larger than K, thus the spectrum of the 
lower membrane alwalys stays below the corresponding 
single membrane spectrum. 



4.4 Real-space fluctuations 

For a calculation of the mean-square displacement (MSD) 
of membrane i, correlations of the Fourier coefficients are 
necessary (cf. equation (|30[) ). Since the system with dis- 
crete tethers is not translationally invariant, all correla- 
tions have to be taken into account, not only those for k 
and k' = — k that comply with the fluctuation spectrum 

<I^T>- 

The height profile of membrane i (i — 1,2) can be 
expressed as 

/i (l) (r) = 2w(r) • c (l) (53) 

by means of 

w(r) T = (l,{cos(q-r)},{-sin(q ■!■)}) . (54) 
Then the MSD can be readily formulated as 
( (r)]V 4 5> r (r)(c<*> cP)w r ,(r) 

r.r' 

r,r' 
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Fig. 8. Typical displacements of the two membranes in model 
F (cf. Fig. [IJ. The scaling functions ur 1 ' (a) and «/ 2 ' (b) are 



displayed. 



is related to the spatial maximum of the RMSD 



for membrane i via equation (|56[l . It is plotted versus the in- 
verse of the dimensionless spring constant KA 2 /m for some 
values of the ratio ki / «2 of the bending rigidities, 
(a): The curve for is situated between its two limiting 
cases. For «i / k,2 —* the case of single tethered membrane with 
the same spring constant is recovered, while for ki/ki — > oo 
the upper membrane in model F behaves like a single tethered 
membrane with half the spring constant. For very stiff teth- 
ers ur 1 ' converges to the value for a single pinned membrane, 
~ 0.124. For weak tethers the discretely tethered system ap- 
proaches the behavior of the corresponding continuous model. 

- (b): The above statements also hold analogously for u> (2) 
of the lower membrane, except that for k\/ki — > oo the lower 
membrane behaves like a single tethered membrane with twice 
the spring constant. — Curves (3), (4), and (7) are identical 
in (a) and (b). The RMSD of a single tethered membrane acts 
as a lower bound for the RMSD of the upper membrane in 
model F and as an upper bound for the RMSD of the lower 
membrane. 



branes tethered to a substrate (cf. Fig. [I]). The two tether 
meshworks are assumed to consist of two quadratic arrays 
of equally strong tethers situated directly on top of each 
other. The RMSD depends on the position in a unit cell 
of the tether meshwork (cf. Fig. 5 in Ref. [TB]). 

For the system at hand, the maximum of the RMSD 
appears in the middle of each unit cell and takes the form 



([^)(r)]- 



knT 



A 



(KA 2 



x wr 



Kl 
K2 



with a scaling function for the upper (i — 1) or lower 
membrane (i = 2). These two functions are presented in 
Fig. [8] They exhibit the following properties in the limit- 
ing cases of rigid (i) or very weak (ii) tethers and in the 
intermediate regime (iii): 

(i) For rigid tethers, K — > oo, the behavior of a single 
pinned membrane is recovered for both membranes. Note 
that this statement only holds if the two involved tether 
meshworks are identical and situated directly on top of 
each other. 

(ii) For very weak tethers, K — > 0, w M approaches the 
result for continuous interactions, i.e., 



1 / Hi 



V& (i) (Kl/«2) 



(57) 



(cf. equations and ([55]) and Fig. 0|. If additionally 
K1/K2 — * (i.e., K2 — > oo if the upper membrane or 
Ki —* if the lower membrane is considered), the behavior 
of a single membrane continously attached to a substrate 
(with the same spring constant) is recovered for both 
membranes. If conversely Ki/k^ oo, both membranes 
act like a single, continuously interacting membrane, but 
with an effective spring constant K e g = K/2 for the up- 
per membrane (since two tether meshworks are arranged 
in series) and an effective spring constant K e g = 2 K for 
the lower membrane (since two tether meshworks — one 
above and one below — are arranged in parallel). 

(iii) For intermediate values of the spring constant a 
cross-over between the two limiting cases (i.e., rigid and 
very weak tethers) occurs around the intersection of the 
limiting curves. In the example of equal bending rigidities, 
ki = K2, the cross-over is situated around KA 2 /k\ ~ 119 
for the upper membrane and around KA 2 /k2 — 53 for the 
lower membrane. Thus, the cross-over roughly remains in 
the vicinity of the cross-over for a single tethered mem- 
brane, KA 2 /k ~ 66 Q3]. 



5 Conclusion and outlook 



where equation (|47j) was used. Herein t(r, i) and t(r',i) 
indicate the components — belonging to the wave vectors 
labelled by r and r' — of the appropriate quadrant of the 
matrix M _1 , i.e., the upper left submatrix for i = 1 or 
the lower right submatrix for i — 2 (cf. subsection 14. ip . 

Aiming at predictions for an experimental implemen- 
tation of model F, we consider the root-mean-square dis- 
placement (RMSD) only for this model, i. e., for two mem- 



A formalism for determining the fluctuation spectrum 
and the root-mean-square displacement (RMSD) of fluid 
membranes coupled by discrete or continuous interactions 
has been presented. Also hybrid models comprising dis- 
crete and continuous interacting potentials are tractable: 
Firstly, if connections between two adjacent membrane 
sheets are built up by both discrete and continuous in- 
teracting potentials, the continuous interaction would lead 
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to additional terms involving the harmonic potential para- 
meter 7. Secondly, membrane sheets with solely continu- 
ous interaction to neighboring objects could be dropped 
in the computation in favor of effective elastics constants 
for the neighboring membrane (s). 

The model of membrane pairs connected by discrete 
tethers and optionally tethered to a substrate exhibit rich 
fluctuation spectra in comparison to the corresponding 
systems containing continuous interacting potentials. In 
certain model configurations a nonmonotonic spectrum 
with two relative extrema (cf. Figs. [5] and [7]) or with one 
relative extrcmum (cf . Fig. [6J (in addition to the prevalent 
relative extremum at the spectral origin) was discovered. 

The calculated RMSD serves as an estimate of the 
typical membrane elongations. This gives usefull hints to 
experimentalists, how the spacings (i.e., the rest-length 
of the tethers) between the membranes or/and the sub- 
strate have to be chosen in order to avoid undesired con- 
tact between them, e. g., if membrane-proteins are present 
which should not adhere to the substrate. Already in the 
continuous interaction model the RMSD of a single teth- 
ered membrane can be increased or decreased by putting a 
second membrane below or above, respectively. Assuming 
equally strong interacting potentials, for identical bend- 
ing rigidities («i = K2) the relative change of the RMSD is 
—0.05 or +0.16 and for an extremely soft lower membrane 
[k\I Ki — > 00) —0.16 or +0.19 for the lower or upper mem- 
brane, respectively (cf. Fig. [4j) . Clearly, an upper interact- 
ing potential stronger than the lower one leads to more 
pronounced changes of the RMSD. While the stiffness pa- 
rameter 7 of two equal continuous interacting potentials 
does not influence the relative change in the RMSD, the 
spring constant of tethers plays a decisive role here (cf. 
Fig. [8|). In short, a two-membrane system could be uti- 
lized to modify typical membrane elongations to a certain 
extend. 

The calculation presented here for two and three inter- 
acting membranes can be generalized to stacks of mem- 
branes 35 . How different tether meshworks (e. g., a mesh- 
work shifted laterally with respect to the lower one) in- 
fluence the fluctuations could be worth analyzing in fu- 
ture studies. Likewise, rather than looking at an effec- 
tively harmonic potential it would be interesting to in- 
clude non-linear interactions describing steric interaction, 
van der Waals interaction, contributions due to hydration, 
and possibly electrostatic interaction [30, 33J . Finally, ac- 
counting for laterally mobile tether attachment sites in 
real biomembranes remains another challenge. 

R.-J. M. thanks R. Richter and G. Majer for helpful 
discussions concerning experimental aspects. 



to three membranes 
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'^/i{^ =i |[V 2 ^(r)] 2 

+ELil[ /l<,,(r " ) -^ +1)(r « ) ] 2 } • ( - :) 

A Fourier ansatz for the three heigth profiles (cf. equa- 
tion ([2])), separation in real and imaginary parts (cf. equa- 
tion ([4])) and transformation to independent variables (cf. 
equation ((6|)) leads to 

H^^LxLy^ [(«i|q| 4 +7i)(< 2 + ^ 
+ Mq| 4 +7i+72)(< )2 + 6; 

2)(< 



(•2) 



( K 3 q 
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1 q 
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+ -L x L y [ 7l C 2 + ( 7l + 72 K 2)2 - 2 7 iCa< 2) ] , 

where the translational mode o'q ' of the lowest membrane 
has been excluded to avoid divergences. 

Pursuing the procedure of subsection 13. 1 1 degrees of 
freedom of two membranes are integrated out, resulting in 
an effective Hamiltonian for the component membrane of 
interest. First, the degrees of freedom of the lowest (i = 3) 
membrane are integrated out, 



a 



(3) 



(.3) 
(•4) 



The effective Hamiltonian for the two upper membranes 
— after a trivial shift of the energy scale — then reads 



— L x L y 
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— L x Ly [ • • - (cf . equation d.2i ) J . 



Now we can focus on the upper membrane (i = 1) and 
integrate out the independent variables belonging to the 
middle membrane (i = 2), 



= / da< 2) 
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what leads to an effective Hamiltonian for the upper mem- 
brane, 



Three membranes connected by continuous 
harmonic potentials 

The Hamiltonian for two continuously interacting mem- 
branes (model B, cf. equation |T])) can easily be extended 
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Consequently, the fluctuation spectrum of the upper 
membrane for non-vanishing wave vectors is (cf. equation 



i^oi 2 : 



kBT/L x Ly 
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For the translational mode (k = 0) this expression would 
diverge, but the exclusion of o,q ' ensures the spectrum to 
be finite at the origin, 



-,(1)12 ■ 
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where the harmonic average of the two spring parameters 
enters (like for two springs in series). 

Now we consider the middle membrane (i = 2) and 
integrate out the independent variables belonging to the 
upper membrane (i — 1) in equation (f75|) . 
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what leads to an effective Hamiltonian for the middle 
membrane, 
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For non-vanishing wave vectors the corresponding 
spectrum of the middle membrane reads 
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The expression for the effective spring parameter 7^ is 
analogous to the effective spring constant of a four spring 
system, where two pairs are connected in series and these 
pairs are connected in parallel. The excluded Oq hinders 
the spectrum at the origin from diverging, 



,(2)|2 • 



k B T/L x L y 
72 



(.15) 



which is independent from the upper membrane. 
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